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Abstract. In this paper we present a scheme for ground state cooling of a flexural 

mode of a nanomechanical beam incorporated in a loop-shaped Cooper-pair box (CPB) 
^ ' circuit. Via the Lorentz force coupling of the beam motion to circulating CPB-circuit 

currents, energy is transferred to the CPB qubit which acts as a dissipative two-level 
i-rt ■ system. The cooling process is driven by a detuned gate-voltage drive acting on the 

$Z| [ CPB. We analyze the cooling force spectrum and present analytical expressions for the 

cooling rate and final occupation number for a wide parameter regime. In particular, 

we find that cooling is optimized in a strong drive regime, and we present the necessary 

conditions for ground-state cooling. 



PACS numbers: 85.25.-j, 85.85. -l-j, 37.10.Rs 
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1. Introduction 

Cooling of a macroscopic mechanical resonator to its quantum mechanical ground state 
has received immense experimental [U [21 El IH 13 E] and theoretical [71 El [9l [TOl [TTl |12l [13] 
attention. The quest for quantum limited control of micro- or nanomechanical oscillators 
is motivated by fundamental interest [HI [15] as well as possible applications in quantum 
information [HI [IH [HI [HI [201 [HI [i] and precision measurement [H [21 [23]. Despite 
these promising and significant experimental and theoretical efforts, to date occupation 
numbers below 1 have not been achieved. 

One can distinguish roughly three routes towards cooling of a nanomechanical 
resonator down to the quantum regime, which have been proposed or partly 
demonstrated: passive optical cooling [31 [H [5l [H [10], active optical feedback cooling 
[H [21 [HI [21] and cooling via coupling to a different heavily damped solid state system, 
such as a superconducting single electron transistor [221 [211 [2Z] , a quantum dot [7] , a 
Cooper-pair box [HI [12] or a flux qubit [13] . 

For passive optical cooling it was shown theoretically in [71 [HI [9] that cooling to 
the quantum mechanical ground state is possible. Experimental work [5] impressively 
demonstrated this cooling method, but did not achieve the ground state due to 
insufficient laser phase noise stabilization. 

For nanomechanical resonators coupled to superconducting circuits less is known 
about ground state cooling, despite recent experimental achievements: It has been shown 
in [221 [23 [2H] that it is possible to integrate a nanoresonator into small superconducting 
circuits and to excite only a single mode of the resonator motion. Also, cooling 
effects have been demonstrated by coupling to single-electron transistors [22] and 
superconducting microwave cavities [29]. On the theoretical side, [HI [H] showed in 
a fully quantum mechanical treatment that it is possible to reach the ground state for 
mechanical frequencies above 100 MHz via capacitive coupling to a Cooper-pair box. For 
other circuit setups a quantum treatment [301 [27] predicts cooling, but final occupation 
numbers beyond the ground state. Cooling effects were also predicted with semi-classical 
approaches in [HI [31] . 

In this paper we perform a full quantum mechanical treatment proving that ground 
state cooling is possible for nanomechanical resonators coupled via Lorentz force to a 
loop-shaped Cooper-pair box qubit. The loop-shaped Cooper-pair box was originally 
introduced to allow for qubit state read-out at an optimal working point [32j, where 
the qubit eigenstates correspond to clockwise and ant i- clockwise currents which can be 
excited and measured. We investigate the conditions for ground-state cooling of the 
fiexural beam motion, using the decay F of the excited Cooper-pair box qubit state 
to dissipate the resonator energy. With carefully chosen qubit-drive detuning, energy 
moves from the resonator motion to the qubit, which acts as an additional reservoir at a 
lower effective temperature. We find that the final occupation number crucially depends 
on the initial thermal occupation A^m of the resonator mode and the dimensionless 
parameter (3 = g'^/iV'-y^, with 7jn the intrinsic dissipation rate of the resonator and 
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g the resonator-qubit coupling strength. Most efficient coohng is achieved for large /3, 
which corresponds to the strong coupling regime of cavity QED. Our analysis shows that 
cooling is optimized in the strong-drive regime for the qubit, and ground state cooling 
is possible even for 10 MHz mechanical frequency. 

Our quantum mechanical model is developed in close analogy to the description of 
laser cooling of trapped ions [331 121] ■ The basic working principle there is that a fast, 
dissipative dynamics of the internal degree of freedom couples weakly to the motional 
degree of freedom and provides an effective dissipative evolution for the latter. This 
is translated here to the solid state context, where the CPB takes over the role of 
internal electronic levels and the mechanical resonator corresponds to the ion motion 
in a harmonic trapping potential, as was done in |71 [11] . Readers not so familiar with 
the quantum optical techniques of master equations and adiabatic elimination based on 
projector operator formalism will find many useful details of calculations in the appendix 
to this paper. 

This paper is organized as follows. In section [2] we present the nanoresonator 
beam and the CPB circuit, and derive the Lorentz force interaction term. Section [3] 
is devoted to a quantum mechanical derivation of the cooling equation, starting from 
a master equation of the full system. In section [H we discuss the cooling rate and the 
final occupation number of the resonator motion in various parameter regimes, before 
we conclude in section [5l 

2. The setup 

In this paper we study a nanomechanical resonator beam integrated into a loop-shaped 
Cooper-pair box (CPB) circuit, with the goal to derive conditions for ground-state 
cooling of a flexural mode of the nanoresonator via the Lorentz force. The CPB circuit 
here acts as a dissipative two-level system (qubit) which removes energy during the 
cooling process. In this section we present a model for the motion of the nanoresonator 
and describe the CPB circuit with the integrated mechanical beam. The system 
Hamiltonian consists of three parts, 

where H^ describes the dynamics of the nanoresonator, ifq represents the driven CPB 
qubit and ifint describes the interaction between the motion of the nanoresonator and the 
CPB qubit. The coupling to the environment will be included using a master equation 
formalism. 

2.1. The nanomechanical resonator 

The nanomechanical resonator [35] considered here is a quasi-lD beam of /x-scale length 
that executes mechanical oscillations. In this setup we consider a so-called doubly 
clamped beam - a beam fixed at both ends - which can undergo longitudinal, torsional 
and flexural oscillations. We are mainly interested in flexural oscillations. Assuming a 
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homogeneous rod, the frequency for the n-th mode of flexural oscillations is given by 
ujn = ^JEMB/{pA)f3'^ [35], where E is Young's modulus, p the density and A is the cross 
section of the beam. Mb = h^w/ 12 is the bending moment for a beam of rectangular 
cross section of width w and height h, oscillating vertically. l3nL is a numerical factor 
which for the lowest mode is f3i ~ 4.73/ L, where L is the length [35]. In the following 
we restrict our analysis to the lowest flexural mode, described by a harmonic oscillator 
Hamiltonian H^, 

H^ = hw^a'^a, (2) 

with a the annihilation operator satisfying bosonic commutation relations [a, a"!"] = 1. 
With the beam integrated into a superconducting circuit, we assume an aluminum 
construction with typical values E = 70 GPa and p = 2.7g/cm^. Choosing a height 
h = 0.3 nm and a length L = l/xm, we get a typical frequency uj^ =10 MHz. Note that 
the frequency does not depend on the width w, but we have to choose w <^ L since we 
treat the nanorod as a quasi- ID beam. 

The dynamics of the nanoresonator mode is further governed by coupling to 
phonons in the electrical circuit and by the residual coupling to all other modes. 
Assuming that the lowest flexural mode is sufficiently well isolated, we model the effect 
of its environment as the interaction with a thermal bath, which can be treated in terms 
of an effective master equation [36] with the Liouvillian £in,d(p); 

CmAip) = 7m(A^m + l)^H(p) + ^rnN^V[a^]{p) (3) 

where 

V [A] (p) = 2ApA^ - A^Ap - pA^A. (4) 

The decay rate 7m = oj^JQ = 0.1 kHz for a typical Q- value of 10^. At cryogenic 
temperatures (T = 15 mK) the steady state of the nanoresonator is a thermal state 
corresponding to an occupation number A^^ = [^wif^ra/ {k-QT)) — 1]^^ ~ 200. 

2.2. The Cooper pair box as a dissipative two level system 

A Cooper Pair box (CPB) [371 ES [39] is a superconducting island connected via one or 
two Josephson junctions ( JJ) to a superconducting reservoir. The dynamics of the CPB 
is governed by the charging energy Ec, given by Eq = (2e)^/(2Cs) with Cs the total 
capacitance of the island, and the Josephson energy Ej, which determines the tunneling 
of Cooper-pairs between the island and the reservoir. The dynamics is restricted to 
controlled Cooper-pair tunneling if the superconducting gap As is larger than the single 
electron charging energy, Ag > -Ec/4, and thermal processes are weak, k^T <C Eq. 



Under these circumstances, it is possible, as we will show in Appendix A to define a 
two-level system (qubit) of the two lowest energy eigenstates of the island, \g) and |e). 
In this paper we consider the split CPB [32], as shown in figured} where two (ideally) 
identical JJs, each with capacitance C and Josephson energy Ej, couple the island to the 
reservoir. Transitions between the two states are induced by a gate voltage microwave 
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Figure 1. Two ways of creating currents in a split CPB circuit. The CPB consists 
of a superconducting island coupled by two Josephson junctions, each with Josephson 
energy Ej and capacitance C, to the superconducting electrode. For a specific gate 
voltage Vg the qubit eigenstates correspond to circulating currents in the circuit. 
Currents are excited (a) by a magnetic flux <i>e applied through the loop or (b) by 
coupling the CPB to a large current-biased Josephson junction. 



drive 6V„ cos uj^j.t, applied via a gate capacitance Cg together with a static bias K 



As we have shown in Appendix A , the truncated circuit Hamiltonian reads, in a frame 
rotating with the microwave drive u^t, 

Ml 



hA 
H^ = — -(T, 



-CTrr 



(5) 



with the detuning A = cudr — '^Ej/h and the drive strength fl = EcSrig/H. H^ is written 
in the eigenbasis at the charge degeneracy point, defined by CgVg,o/(2e) = 1/2 and 
drig = CgdVg/ {2e) = 0. At this particular bias point the system is less sensitive to 
charge fiuctuations [321 HO], which would transform into fiuctuations of the detuning A 
and in principle could make the cooling procedure less efficient. For a more realistic 
setup, one has to assume an asymmetry of a few percent; we will discuss this issue in 



Appendix A 



As already mentioned in the introduction, the loop-shape of the split CPB allows 
for circulating currents in the circuit. In fact, the energy eigenstates at the charge 
degeneracy point correspond to clockwise and anticlockwise loop currents [1^. In the 
CPB circuit as described so far there is no asymmetry in the loop allowing the fiow of 
a net current in either direction and therefore we have to add elements which break the 
symmetry. One possibility is to thread the circuit loop with an external magnetic fiux 
$e, as shown in figure [l](a). The magnetic field identifies a direction in the loop, and 
consequently one current direction is energetically favorable, giving rise to net currents in 
the circuit. The other option is to integrate a large current-biased Josephson junction 
with Josephson energy Eji, into the loop, as shown in figure [T](b). As we show in 



Appendix A[ the bias current Jb determines the phase across the large JJ, and thus acts 
in a similar way as the magnetic fiux $e. For the purpose of this paper, both ways to 
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excite the loop currents are equivalent, and give rise to loop currents ±Jq , 

with 00 defined by $c = {h/{2e))(j)o and sin</)o = (V(2e))/b/-E'j,b5 respectively. It is 
obvious from the expression ([6]) that in the absence of magnetix flux or bias current, 
00 = and the loop currents vanish. Apart from the excitation of loop currents, the 
energy splitting is affected and thus the detuning A is modifled, 

A ^ u;dr r- ^os ( Y ) • ( ' ) 

As a last remark, we would like to point out that this derivation holds also for Eq ~ ii^j, 
where the qubit is even less sensitive to charge fluctuations |32l HO]. Due to the 
anharmonicity of the energy spectrum, it is possible to truncate the Hilbert space to 
the two lowest energy eigenstates at the charge degeneracy point 5n^ = 0, and it can 
be shown that these states are eigenstates also of the truncated loop-current operator 
[HlllJ]. However, more charge states will be involved in the qubit eigenstates, and thus 
the parameters A and Q are modifled. 

So far we have only described the coherent dynamics of the CPB. Due to its 
coupling to the surrounding circuit leads and uncontrolled degrees of freedom within 
the circuit 00], the qubit is subject to dissipation. The dissipative dynamics of the 
qubit is determined by the depolarization time l/(2r), which describes the relaxation 
to the ground state, and the dephasing time l/(r + Fd) [40j. These two processes are 
included phenomenologically [39] into the qubit master equation, 

C^M = TV[a_]ip) + ^V[a,]ip) (8) 

with V[a_]{p) and V[az]{p) given by equation (jl]). For an energy splitting Ej/h ~ 20 
GHz [32] at a cryostat temperature T of 15 mK, the thermal occupation number of 
the excited state is negligibly small, A'q = [exp^Ej/^ksT)) — 1]^^ ~ 0, and therefore 
the heating term NqT>[a^]{p) can be neglected. Both F and Fj can vary dramatically 
between different samples and working points [15| [32] . Since the relaxation is crucial 
for the cooling procedure, we assume a typical value F ~ 1 MHz [32]. The derivations 



presented in section |3l Appendix B and Appendix C are valid for a flnite Fd, but for 



transparency we restrict the discussion in section H] to the limit Fd/F — > 0. 

2. 3. Lorentz force interaction 

Now we include the nanoresonator beam into the loop of the CPB circuit as shown in 
figure[2]^a). As already mentioned earlier in this section, it is possible to excite circulating 
currents in the circuit, and we assume a construction where the currents are also driven 
through the resonator beam. A uniform magnetic field Bq applied perpendicular to 
the long axis of the beam, as shown in figure [2t^a), together with the circulating qubit 
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Figure 2. (a) Inclusion of the nanoresonator into the CPB circuit. The coupUng 
between the two systems wih be provided by the apphcation of a magnetic field along 
the y-axis in plane with the CPB, such that the current flowing in the circuit along 
the x-direction will couple to the motion of the oscillator via the Lorentz force, which 
acts in the z-direction as sketched in (b). 



currents /q ([6]) creates a Lorentz force Fl which acts on the beam in the perpendicular 
direction z [351 El HH] as shown in figure El^b), 

Fl = bJl. (9) 

Since the sign of the Lorentz force is qubit-state dependent through the qubit operator 
/ (1A.12I) . it can act as a coohng force on the displacement of the beam in the z-direction 
through the interaction term ifint = Fi^z. Writing the displacement of the beam as 
z = ao{a + a"'"), with oq = y/h/{2muj^) the amplitude of the zero-point motion, the 
interaction Hamiltonian gets the form, 

^int = hga^i^a + a''), 

with the coupling constant g given by, 
BoIqLao 



(10) 



9 



h 



:ii) 



Note that this interaction term differs from the ubiquitous Jaynes-Cummings interaction 
|46j , and thus the requirements for ground-state cooling will differ from for example the 
case of ion cooling [33]. In this paper we focus on the case of constant g. Assuming 
ao ~ 10~^^ m, Bq f» 10 mT and Jq ~ 10 nA, we achieve a typical coupling constant 
g ^ 100 kHz. In principle it would also be possible to drive either the coupling magnetic 
field Bq or the bias current /b, but this is beyond the scope of this paper. Note that the 
coupling field Bq does not create an additional flux in the circuit loop. However if we 
choose to use the setup in flgure[T|^a), the applied flux $e in the i-direction would create 
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a Lorentz force in the perpendicular direction y and could in principle lead to heating 
of the oscillation modes in this direction. However this coupling is very weak, since the 
magnetic field required to obtain $e is at least one order of magnitude lower than the 
coupling magnetic field Bq. In addition, assuming a rectangular beam with w ^ h, the 
oscillation frequency for the motion in the y-direction would be large compared to cUm- 

3. Cooling 

In this section we derive an effective master equation for the relevant mode of the 
resonator motion, and the corresponding cooling equation. Following the line of |34j . 
our main assumption is g <^r, which means that the qubit goes to equilibrium on a time 
scale which is much shorter than the time scale for the interaction between the qubit 
and the resonator. Thus the state of the qubit is hardly affected by the interaction, 
and the qubit can be adiabatically eliminated. The second important assumption is 
7m ^m ^ T- Due to the different time scales for the dissipation of the resonator via 
the qubit and via the coupling to its own environment, the two dissipation channels 
act separately. We define the small parameters ei = g/T and 62 = 7mA^m/r, and 
thus treat the coupling between the qubit and the nanoresonator and also the intrinsic 
dissipation of the resonator as perturbations. For the perturbation expansion of the 
master equation, we write it on the following form, 

p = C^ip) + CM + ^lip) + Mp) (12) 

where the zeroth order terms describe the non-dissipative dynamics of the resonator, 

C^{p) = -^[H^,p] (13) 

and the dynamics of the qubit as described in section [2l 

cM = -1[h^,p] + ^^Ap)- (14) 

The first order term in ei describes the coupling between the qubit and the resonator, 

Ci{p) = -^[H,^up] (15) 

whereas the term of first order in €2 gives the coupling of the resonator to its environment 
as described in section [21 

Ciip) = >Cm,d(p)- (16) 



We shift all details of calculations to Appendix B where we derive the reduced master 



equation for the motion of the resonator, and state here directly the result, 

p^ = -iLJ^[a^a,p^] + A_V[a]{p^) + A+P[a^](pm) (17) 

with pm = trq{p} given in a frame where the resonator motion is coherently shifted as 
defined in ( ]B.1|) . with the shift a given by (1B.9|) . The effective frequency is given by 

OJm = ^m + g^lm{S{u^) + S{-UJm)}, (18) 
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where the frequency shift is an order ~ g'^/{uj^^r) ~ 10"'^ smaller than the bare resonator 
frequency cj^- The cooling and heating rates are given by A^ and A^ respectively, 

A_ = 2 [g^Re{S{uj^)} + 7m(iVm + 1)] , A+ = 2 [g^Re {Si-cu^)} + 7miVm] • (19) 

Here we defined the spectral function 



S{uj)= I dTe'^^{Sa,{T)5a,),,, (20) 



where Sa^ = Oz — {c''z)ss is an operator describing fluctuations of a^ about its mean 
steady state value (cr2)ss- The coupling to the CPB qubit has essentially two effects on 
the mechanical resonator, a negligibly small frequency shift and an additional dissipative 
channel, which can eventually exceed the natural dissipation due coupling to a phononic 
heat bath. If paramters are chosen in the right way, this can give rise to strong cooling. 
Note however that the master equation for the resonator motion (TT7j) is correct only up 
to 0(e2(,)). 

From the master equation we derive the corresponding equation for the expectation 
value of the motional occupation number (n), 

j^{n) = -{A_~A^){n)+A^. (21) 

The cooling equation ( 12T1) gives us information about the total cooling rate W, 

W={A^-A+), (22) 

and, provided that we are in the cooling regime {W > 0), the flnal occupation number 
is given by 

(n), = ^. (23) 

Equations ( l22l) and [23] are the main results of this paper. 

4. Discussion 

In the following, we will analyze the cooling rate fl22|) and the flnal occupation number 
( !23l) . and deflne parameter regimes where cooling to the ground state is possible. In 



Appendix C we derive an expression for the spectrum S{uj) ( JC.71) which provides general 



expressions for the cooling rate W and for (n)f in terms of the system parameters 
A, i7, r,7m,co'm, A^m, gf. Wc use the full expression to generate plots for {n)f and W 
for a realistic set of parameters which were motivated in section [2l Furthermore we 
give analytic expressions in the limit of weak qubit decay F <^ A, where we introduce 
the qubit energy splitting in the rotating frame A = i/fi^ + A^ - see the single qubit 
Hamiltonian ([5]). 

The cooling rate consists of two parts, as the coefficients A± (IT9l) show, where one 
part is due to the intrinsic decay 7ni of the resonator and the other one is due to the 
coupling g to the dissipative qubit. The qubit serves as an additional bath for the 
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Figure 3. Real part of the force spectrum (in units l/F) for the parameters: Wm = 10 
MHz, r = 1 MHz, A = 10 MHz, Vl = 4MHz (dashed, red), Vl = 8.5 MHz (sohd, green) 
and Vt — IQ MHz (dotted, blue). The dashed curve shows the absence of the heating 
peak for weak drive, but still slow cooling due to the small amplitude of the cooling 
peak. The solid line shows the optimal cooling, where the cooling peak is high and the 
heating peak still small. For even stronger drive, represented by the dotted line, the 
heating peak has increased so much that cooling and heating balance each other and 
the net cooling effect is lost. 



resonator, with an effective decay rate determined by the couphng strength g and the 
real part of the force spectrum S{(jj) (l20l) . Thus it is the presence of the Lorentz force 
which makes it possible to cool the resonator motion below its thermal occupation. In 
the limit of weak qubit decay rate, F ^ A, the real part of the spectrum exhibits well- 



resolved peaks at ±A as we have shown in Appendix C To lowest order in F/A the 
spectrum in the vicinity of the peaks can be written. 



Re{^(cj~±A)} 



r+a± 



n^ 



a± 



T 



A ^2 



(24) 



F^ + (cutA)2' "^ 2A2 ' A2A2 + fi2^ 

with the effective decay rate F+ presented in (JC.lOp . Since only the amplitude of the 
real part of the spectrum at the resonator frequency. Re {S'(±co'm)}, enters the cooling 
equation (|2T|) . it is intuitively clear that in this regime cooling is significant only for 
A = cUm- Note that Re{S'(ci;m)} is responsible for cooling, whereas Re {5'(— tUm)} 
contributes to heating. When inserting equation (p^ into the cooling rate (122|) . we 
observe that a net cooling effect requires a+ > a_, i.e. negative detuning A < 0, which 
will be assumed in the following. We plot the real part of the spectrum as a function 
of the dimensionless frequency uj/T in figure [3] for A = ujja and different values for Vt. 
The dashed curve for which we chose Vt = 0.4 tUm shows the absence of a heating peak 
for weak drive, but still slow cooling due to the small amplitude of the cooling peak. 
The solid line {Q = O.SStUm) shows the optimal cooling, where the cooling peak is high 
and the heating peak still small. For even stronger drive {Q = cUm), represented by the 
dotted line, the heating peak has increased so much that cooling and heating balance 
each other and the net cooling effect is lost. 

The total cooling rate as a function of detuning A and drive strength Q is shown in 
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figure 111(a). We observe that cooling is obtained for a negative detuning (A < 0), as we 




Figure 4. Contour plot of (a) the cooling rate (in units of kHz) and (b) the final 
occupation number as a function of the detuning A and the drive strength O for the 
following choice of parameters: A^m = 195 , 7,11 — 17 Hz, F = 1 MHz, Wm = 10 MHz, 
g = 0.1 MHz. The dotted line shows the special value 57^ + A'^ = Wj„, where cooling is 
optimal. 

also have seen in the calculation of the spectrum. This means that the energy provided 
to the qubit is less than its excitation energy, and thus energy is taken from the motion 
of the resonator to fully excite the qubit. Consequently the resonator will be cooled. In 
addition we notice that the highest cooling rate is determined by a resonance condition 
A = ujja as required by energy conservation. This is clearly visible from the dotted line 
in the A — fi-plane of figure|l](a). A similar resonance condition has been observed for an 
LC-oscillator inductively coupled to a flux qubit in [17] . This result is also obvious from 
the analytical derivation of the spectrum ( 12^ performed in Appendix C Thus there is 



a trade-off between the optimal detuning and the optimal drive strength. Furthermore 
it is clear from figure Hl^a) that the region of optimal cooling covers only a small section 
of the dotted line A = u^. This motivates us to introduce the angle (p defined by 

sin(^ = fi/A. (25) 

We find that the highest cooling rate {W > 8 kHz) appears in a region around 
simp ^ 0.85, which is confirmed analytically below. This means that fl has to be 
remarkably different from zero in order to get a high cooling rate, i.e. we have to drive 
the system strongly. We also note that we cannot cool the resonator below its thermal 
occupation for A = or for ^2 = 0. The latter is easily understood since we need 
to drive the qubit in order to bring it into the excited state from where it can decay. 
For A = the qubit is driven at resonance, and thus will not take energy from the 
resonator motion. This can also be understood from the the dotted spectrum in figure 
[31 where the positive and negative frequency parts are equally high; the spectrum ([JT 
gives Re{5(±cjm)} = l/(3r) for A = 0. 
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Figure 5. (a) The dimensionless function /((/j) plotted as function of sunp = ri/wj, 



with 9? + A^ 



The maximum is obtained for svmp = O.SSwm- (b) The final 



occupation number as a function of the drive strength sin Lp = 57 /wm for a specific set 
of parameters: cj,n = 10 MHz: g = 0.1 MHz, 7^ = 17 Hz, T = 1 MHz, iV^ = 195. 
The minimum is in this specific case achieved for sin<y9 = 0.80, with the minimum 
occupation number 0.8. 



Now we study the cooling rate in more detail analytically, in the limit of weak qubit 
decay and optimal detuning (F <^ A = uj^. Inserting the spectral amplitudes (12^ in 
the cooling equation (I2T1) . we derive the total cooling rate W , 

W = A_-A+ = 2^^[l + [3f{ip)], (26) 



with the dimensionless parameter /5, 

/3 



9' 



7n 



(27) 



4sin^ V^V 1 — sin^ 



<^ 



/ [4 — sin'' ip\ which is plotted as a 
\p) describes the behaviour along the 



and the function /(v?) 

function of sinyj in figure [5]^a). The function f\ 

dotted line in figure ll](a). The analytical expression for the cooling rate ( 126|) tells us 

that fast cooling is obtained for a large value of the parameter (3 and for a drive strength 

Vt which maximizes the function /(</?): /max = /(^o.w) with f2o,w = 0.85 tUm- Following 

the dotted line in figure Hl^a) and the plot of the function f{ip) in figure [5](a), we see 

that the cooling rate increases very slowly for weak drive, 



PF ^27„ 



1+/3 



fi< lAI 



(28) 



up to the maximum point f2o,W) from which it quickly drops to its minimum value 

W = 27m for very strong drive, 

'4|A| 



ly ~27„ 



l + (3 



sn 



n 



UJ„ 



(|A|-0). 



(29) 



Figure 111(b) shows the final occupation number (n)f as a function of the detuning A 
and the drive strength Q for an initial thermal occupation of A^m = 195. The minimum 
final occupation number {n)f < 0.85 for this specific choice of parameters is also obtained 
for the detuning A = uj^ which follows from energy conservation represented by the 
dotted line. The optimal drive strength is in the region where siny9 ^ 0.8. Although the 
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final occupation number is not exactly proportional to the inverse of the cooling rate, 
for the parameters we have chosen here the minimum of {n)^ almost coincides with the 
maximum of the cooling rate, which is promising for efficient cooling procedures. 

Now we study the final occupation number analytically in the limit F <C A = ti^in- 
By combining the spectrum (l24l) with the cooling equation (12T1) we derive the final 
occupation number (n)f, 

A+ N^w+{^)+(3a_ 2 + sinV .^^x 

(n)f = — — = -, wAiD) = . 30) 

From the expression fl30p it is clear that for zero driving {simf = 0), the final occupation 
is equal to the initial (thermal) population, {n){ = N^, and thus there is neither cooling 
nor heating. For zero detuning [simp = 1) the driven process heats up the system 
and finds a steady state at {n){ = N^ + (3/3. Figure \5Kh) shows {n)f as a function of 
sin (fi = Q/A for a specific set of parameters and thus gives a feeling of the qualitative 
behaviour. Following figure [H](b) or the dotted line in figure |H^b), we see that the final 
occupation number decreases slowly for weak drive, 

and finds a minimum close to the maximum for the cooling rate (here at r2o,Nf = 0.80 Uj^) 
after which it quickly increases again, 

Wf- 77777\' 1^1 ^0. (32) 



1 + K^ 



In contrast to the optimum drive strength f2o,w for the cooling rate W, the drive strength 
Qo,m which minimizes the final occupation number depends on the parameters A^m and 

Assuming optimized drive strength n^Nf, the thermal population N^ and P are 
the only two parameters determining the final occupation number. These parameters 
represent two competing processes; the intrinsic dissipation of the resonator represented 
by N^ and 7m, and the driven cooling of strength g'^/T. For the latter process to be 
dominant, it is required that (3 is significantly larger than N^, while respecting the 
basic assumptions g/T ^ 1 and 'jmNm <^ F of the derivation. In figure [6] we study 
the behaviour of the final occupation number beyond the limit of resolved sidebands. 
Figure [6](a) shows how the minimized final occupation number changes with the thermal 
resonator occupation N^ and the universal parameter /3 = g'^/{'y^T). We clearly see 
that the larger temperature of the resonator (i.e. larger N^), the larger (3 is required to 
reach a low final occupation number. A larger f3 is obtained either hj stronger resonator- 
qubit coupling g or by weaker decay rates F and 7m. For N^ = 195, corresponding to 
Ura = 10 MHz at T =15 mK, ground state cooling is possible provided that P > 500. 
In figure [6](b) we study the final occupation number as a function of the quality factor 
Q = oJrnllm of the resonator and the decay rate F of the qubit, for a specific thermal 
occupation number N^ = 195 and fixed coupling g = 0.1 MHz. A lower occupation 
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Figure 6. Contour plot of the final occupation number (a) as a function of the 
universal parameter j3 — 3^/(7mr) and the initial occupation number of the nanorod 
A^m and (b) as a function of the quality factor Q of the nanorod and the decay rate F 
of the qubit for the following choice of parameters: N^i ~ 195, g = 0.1 MHz. In both 
plots each point has been optimized for the detuning A and the drive strength Q . 



number is achieved by increasing the quahty factor of the resonator, or by decreasing 
the decay rate of the qubit. For F = 1 MHz ground state coohng can be obtained for 
Q > 5 ■ 10^, which corresponds to 7m < 20 Hz for tUm = 10 MHz. This corresponds to 
the resolved sideband hmit and the parameters we used so far. 

Remember that the resuhs are presented in a frame where the resonator motion 
is shifted by a coherent shift a ( IB. 91) . In the original frame we would have to add the 
coherent part ~ |ap to (n)f. This shift can be removed by adiabatically switching off 
the coupling magnetic field Bo or the circulating currents Jq. 

The results presented in this discussion are valid in the limit of negligible pure 
dephasing, Fd/F -^ 0. For finite Fd one has to consider the more general analytic 
expressions in Appendix C for the cooling rate W (JC.121) and the final occupation 



number (n)f (10.130 . which are valid for F,Fd <^ A. From these expression we draw 
the conclusion that our results are still valid for weak dephasing rate Fd -C F, whereas 
a large Fd would change the conditions for ground-state cooling considerably. 

As a last remark, we point out the similarities and differences of the cooling process 
discussed in this paper compared to the laser cooling of a trapped ion [33] . The dynamics 
of the two-level system in this setup (the CPB qubit) is formally equivalent (in the limit 
Td/T -^ 0) to the damped and driven two level system modeling the ion, with its 
dynamics completely described by the Bloch equation (IC.ip . The position of the peaks 
in the spectrum is fully determined by the eigenvalues of the matrix A ( 1C.2I) . This 
explains why we obtain optimal cooling in the same regime, where cooling for a trapped 
ion is optimal, namely in the resolved sideband regime. The main difference between 
the two systems is the cooling force, which in this setup gives an interaction of the 
form ~ a^(a + a^), in contrast to the Jaynes-Cummings term ~ (c+a + a^a'') which 
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appears for ion cooling. This difference manifests itself in the dependence on A and VL 
of the amplitude and the effective width of the spectral peaks, and consequently affects 
the optimal drive strength. As already mentioned earlier, there are contributions to 
the cooling rate both from the intrinsic dissipation 7in of the resonator and through 
the coupling to the two-level system, in contrast to ion cooling where only the latter is 
present since the motion of the ion is perfectly isolated. Finally we should mention that 
whereas the laser both drives the internal degrees of freedom of the ion and its coupling 
to the motion, in this setup the qubit coupling to the motion is independent from the 
microwave signal which drives the qubit. 

5. Conclusion 

In this paper we have derived the conditions for ground-state cooling of a nanoresonator 
beam which is integrated in a loop-shaped Cooper-pair box (CPB) circuit. The CPB 
qubit acts as leaky two-level system through which the energy of the resonator motion 
dissipates. The cooling force in question, the qubit-state dependent Lorentz force, acts 
on the beam due to circulating currents in the circuit in combination with an applied 
magnetic field. We have derived a general expression for the force spectrum, which 
determines the cooling rate and the final occupation number, and evaluated the result 
in the resolved sideband regime where the qubit decay F is small compared to the 
resonator frequency oj^a- 

In the resolved sideband regime, the spectrum exhibits narrow peaks and cooling 
is optimized for the detuning yAM-TF = uoja- Assuming this condition, we further 
investigated the cooling rate and the final occupation number as a function of the drive 
strength fi, and found that optimal cooling requires a strongly driven qubit; the cooling 
rate is maximal for Q, = 0.85 c<Jm independently of the other parameters. The final 
occupation number is minimized in the region Q ~ 0.8 uj^ for our parameter example. 
When optimizing the final occupation number with respect to the drive strength, we 
found that basically two parameters determine the possibility to cool the resonator 
mode to the ground state, namely the initial thermal occupation number N^ of the 
resonator, and the parameter j3 = g"^ /{T'-)^. Thus given an initial occupation number 
of the mechanical mode, the possibility to reach the ground state is determined by the 
relation between the coupling strength g and the decay rates F and 7m, requiring strong 
enough coupling in comparison to the decay rates. Generally cooling requires (3 > N^, 
and in particular we could show that for N^^^ = 195, corresponding to uj^ = 10 MHz at 
T = 15 mK, ground state cooling requires (3 > 500. 

Let us also say a few words about the limiting cases. For = 0, the qubit is 
not driven and provides neither cooling nor heating. The cooling rate is given by the 
dissipation rate of the resonator, 27^, and thus the mechanical mode stays at thermal 
occupation (n)f = Nj^. For zero detuning, A = 0, the qubit provides both cooling and 
heating of the resonator motion with equal amplitude, and the resonator is heated up 
to the final occupation number N^ + f3/3. The total cooling rate is still given by 27m. 
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From the form of the parameter j3 it is clear that for a specific couphng strength 
g, coohng is improved by a larger Q of the resonator, which basically reflects the fact 
that the resonator environment needs longer time to heat up the mechanical mode to 
its thermal state after it has been cooled. 

Cooling is also improved by a weaker qubit decay F, which is not very intuitive 
since the cooling procedure relies on the dissipation of the excited qubit state. However, 
from the point of view of the resonator, the decay of the qubit state shows up as an 
effective decay rate g'^ /T in for example the heating and cooling rates A±, and therefore 
the cooling efficiency decreases with larger F. 
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Appendix A. Derivation of the Hamiltonian for the Cooper pair box 

In the following we derive the qubit Hamiltonian ([5]) for the CPB circuit shown in figure 
[D^b). Let us start from the circuit Lagrangian L(0j, 0j) [H], which is presented in terms 
of the phases 0i,02,0 across the respective Josephson junctions (JJs) of the circuit, 
using the voltage-phase relation V = (^/(2e))0, 

L = LscB + Ln, (A.l) 

where the first term describes the driven Cooper pair box. 



^SCB - I - 



^(0? + 0^) + ^ (1^.-0" 



+ Ej(cOS0i + COS02) (A.2) 



and the second term describes the large current-biased Josephson junction (JJ), 

-J ^02 + Ej,bCos0+-/b0. (A.3) 

In the following we assume the phase regime E'j^b ^ -E'cb = (2e)^/(2Cb) for the large 
JJ, in which case the dynamics of the phase is similar to that of a particle in a tilted 
washboard potential, with the tilt determined by the applied bias current Jb. Provided 
that the tilt is not too large, /b < /c where Ic = (2e/fi)£'j,b is the critical current of the 
JJ, the potential contains well-defined wells in which the phase is trapped. Further we 
assume that Cb — >■ 00, corresponding to a particle with infinite mass, in which case the 
phase stays at the minimum point 0o determined by the bias current via the relation 
sin 00 = Ih/h^ and the large JJ can be treated as a classical device without dynamics. 
Here we also assumed that the currents through the large JJ due to the CPB dynamics 
are small, which is justified if Ej/Ej^i, ^ 1. The contribution from the large JJ to the 
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Lagrangian is thus only a constant term -Ljj(0 = 0o)- On the other hand, assuming 
neghgible self-inductance in the loop, the phase 00 determines the phase across both 
CPB junctions through the flux quantization relation [49], 

01 + 02 + 00 = 0. (A.4) 

A similar effect is achieved by applying a magnetic flux $e as shown in figure [It^b). In 
the following it is understood that 0o either is determined by a bias current or by an 
applied magnetic flux with 0o = (2e//i)$e- We are left with only one degree of freedom 



(A.5) 



+ 2Ejcos(^) cos0q. (A.6) 

In order to derive the circuit Hamiltonian, we calculate the charge number on the island 
riq, which is the conjugate variable to the phase 0q, 

and obtain the circuit Hamiltonian H{nq, 0q) = 'n.q0q — L [50], 

H = Ec K - rig)' - 2Ej cos f y j cos0q (A.8) 

where the charging energy Eq = (2e)^/(2Cs) contains the total island capacitance 
Cs = 2C + Cg, and n^ = CgVg/{2e) is the dimensionless charge induced on the island by 
the gate voltage. Note here that the two-junction circuit effectively behaves as a single 
junction CPB with tunable Josephson energy 2Ej cos(0o/2). Now we can also derive the 
current flowing in the circuit. The loop current is determined by the time derivative of 
the charge q through both JJs, which is conjugate to the phase — (0i + 02) = 0o across 
the CPB. Using Hamilton's equations of motion [50], we get. 



in the system. 






^ 02-01 






0q 2 ' 




and the circuit Lagrangian gets the form. 




\2eJ 


[o>;+ ; 


(2e ■ V 

^\/g+0q 

^h J 



2edH 

h 900 

with the current Jq given by ([6]). Further we quantize the circuit Hamiltonian 



n ocpo 



H = Ec (nq -rigf - 2Ej cos ( y 1 cos 0q (A. 10) 

by imposing the canonical commutation relation [0q, nq] = i. In the absence of tunneling 
{Ej — > 0), the charge states \2en) representing integer number of Cooper pairs on the 
island form an eigenbasis for the Hamiltonian flA.lOl) . For dominating charging energy, 
Eq ^ Ej, the two lowest energy eigenstates are superpositions of the two charge 
states \2en) and |2e(n -|- 1)) [38]. The superposition states are controlled by tuning 
the dc-part of the gate voltage V^^q. In particular, we can bias the island at the charge 
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degeneracy point n^ = 1/2 where \2en) and \2e{n + 1)) are degenerate and the energy 
levels are split up by the weak tunneling. The eigenstates are given by the symmetric 
and antisymmetric superposition states \e/g) = (|2en) ± |2e(n + 1))) /-\/2- Since the 
spectrum is highly asymmetric, it is allowed to consider only these two levels and treat 
the system as a qubit j38], provided that the temperature is sufficiently low. At the 
charge degeneracy point the truncated Hamiltonian for the qubit reads, 

H=ja, + m{t)a^, (A.ll) 

written in terms of usual Pauli matrices in the basis {|e), \g)}. Here we assume 
the energy splitting E = 2Ej cos(0o/2) to be tuned adiabatically, whereas fi(t) = 
(Ec/h) (ng(t) — 1/2) is driven by the gate voltage; n^it) = l/2+6ng coscudrt. Moving to a 
frame rotating with the microwave frequency Udr, we end up with the qubit Hamiltonian 

©• 

During quantization and truncation, the expression for the loop current (jA.QO 

transforms into a current operator which is diagonal in the qubit eigenbasis {|e), \g)} 
[32] 

/ = V.. (A.12) 

This holds also in the regime Eq ~ Ej pT| WI\ . 

So far we have assumed the ideal case of a symmetric Cooper-pair box circuit, i.e. 
one where the two JJs are equally large. In reality the deviation 6E] = {Eji — -Ej,2) is 
usually a few percent, 6E] <^ E;^. We expand the qubit Hamiltonian H (lA.llI) to first 
order in SEj/Ej, 

— o-^ + mX[t)a^ + — 

with 8E = 6Ejsm{(f)Q/2). Further we investigate the effect of the asymmetry on the 
qubit-resonator interaction (TTOl) . In addition to the excited loop current / (]A.12p . the 
asymmetry creates an additional qubit-resonator interaction term via the current /a 

4 = hcTy with 4 = -— cos y , (A. 14) 

which is non-zero even for zero bias current (0o = 0). However, in the rotating frame of 
equations (jSD and (ITU]) , the asymmetry terms are rotating fast, 

(Tj, ^ 2 [e-^^'^'-V_ - e^"*V+] , (A. 15) 

and thus within the rotating wave approximation, a small asymmetry 6Ej/Ej <^ 1 does 
not affect the qubit dynamics or the qubit-resonator interaction. 

Appendix B. Master Equation 

In this section we derive an effective master equation for the motion of the resonator 
starting from the master equation for the full system given in equations flT2l) - flT6l) . First 
we perform an adiabatic elimination of the qubit degrees of freedom assuming T ^ g 



H ^ —ij^ + hQ(t)ax + ^:;-(Ty, (A.13) 
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such that the qubit goes to a steady state on a much shorter time scale than the time 
scale for the interaction of the resonator with the qubit. Thus we can assume that the 
density operator p never deviates much from the factorized form pss ® trqjp} where the 
qubit is in its steady state pss, and project the master equation onto this form. To 
remove the action of the coupling term on the projected state, we begin by defining a 
displacement operator D which performs a coherent shift of the oscillator displacement, 

DaD^ = a + a. (B.l) 

D is a unitary operator which acts only in the resonator space. We define the density 
operator p in the shifted frame, 

p = DpD\ (B.2) 

and similarly applies the displacement operator to the master equation, 

p = Ctp + D{C^p + C,p + C2p)D^ (B.3) 

with 

DC^pD^ = C^p — iij„i[aa'^ + a* a, p], 

DCipD^ =Cip~ig{a + a*)[a,,p], (B.4) 

DC2pD^ = C2P + 'ym[oi*a — aa"^ ,p]. (B.5) 

The shift acquired in DCiD^ acts only on the qubit, and we include it in the qubit 
dynamics by defining a new Liouvillian C^, 

^oP = ^oP - Wia + a*)[(^z, p]- (B.6) 

Then we can write the master equation in the shifted frame as, 

p = (£" + £[J + £1 + C2)p + [(7m - iuJm)oi*a - (7m + iuj^)aa\ p]. (B.7) 
The next step is to project the master equation onto the factorized form, 

PP = Pss ® Pm, Pm = trq {p} (B.8) 

using the projection operator P defined by PC^p = and the projection on the 
orthogonal subspace Q = 1 — P. The shift a is determined by requiring that 
the projection on P of the last term in (1B.7P cancels the projected coupling term 
PCiPp = -igpss ® {(Tz)[a + a^Pm], 

a = - -^^^^. (B.9) 

UJm - «7m 

The master equation in the P-space now reads, 

Pp = PC^Pp + PCiQp + PC2PP, (B. 10) 

where we made use of [P, £™] = and [P, C2] = 0. Remember that our goal is to derive 
a master equation for pm which is correct up to the second order in the small parameters 
ei(2), as we outlined in section [3l The procedure is to derive a closed equation for Pp. 
We first have to replace Qp in f IB.lOp by an expression which only depends on Pp and 
is correct up to first order in the small parameters. This is sufficient as Qp in (IB. 101) is 
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multiplied by £i, which is already of order 0(ei). Using QC^Pp = 0, QC^Pp = and 
QC2PP = we write the master equation in the Q-space 

Qp = Q [C'^ + 4 + £2) QP + QCiPp - igQ[{a, - (a,)) (a + at), Qp]. (B.ll) 

As the term Qp is of order 0(ei), the last term in (]B.11|) is of order O(e^) and the term 
QC2QP is of order 0(eie2). Thus we can neglect these terms and solve the equation 
formally, 

t 

Qp{t) = e«(^S^+^o)(*-*i)Qp(ti) + I dTe'^'^''^^^''^^^'~^^QC^Pp{T). (B.12) 

This equation is exact up to 0(ei(2)). Next we solve the equation for Pp ( IB. 101) in order 
to insert it into the term QCiPp in equation ( lB.12p . To be consistent we have to keep 
terms only up to order 0(e5/2)); 

pp{t) = e^^S^(*-*»)Pp(to). (B.13) 

Being interested in timescales larger than l/F for which Cq decays, the first term in 
( JB.121) will not contribute to the master equation. The integrand in the second term of 
flB.12p decays on a time scale l/F, allowing a Markov approximation where we extend 
the lower limit of integration to ti -^ —00 and replace p(to) by p{t). Inserting the 

solution for Pp(t) (IB. 131) into the solution for Qp(t) ( IB. 121) we obtain, 

00 

Qp{t)= /"(ireO(^S^+^o)-Q£^e-^^S^^Pp(t). (B.14) 


Further inserting this solution for Qp{t) flB.14p into fIB.lOl) . we finally obtain a closed 



equation for Pp. Using the definition for P ( JB.SP and tracing over the qubit degree of 
freedom, the master equation for the resonator motion p^ is expressed as, 

p^ = {Q + C2)Pm + trq < /" drCiQe^''o'+^o)rCiPss ® e-^o"^p^ I . (B.15) 

There is still some work to do to evaluate the second term of the master equation for 
the resonator (IB. 151) . We define the operator, 

ScTz = (Tz- {(rz)ss, (B.16) 

describing fluctuations of cr^ about its mean steady state value {(yz)ss- Further we use 
the following relations: trq < Sa^e^o'^dazPss \ = {So-z{T)6az)ss, trq < Sa^e^o'^pssSaz > = 
{SazSaz{T))ss and e^o^'^a = a{—T). After some algebra we obtain a master equation 
for pm where the qubit degree of freedom is completely traced out. 



{C^ + £.2)pm -9^ dT{6az{T)6az)ss[a + a\ [a{-T) + a^ 



-T),Pn 



/ / rfr([(5a,(r),(5aj)ss[a + at,p„,(a(-r) +at(-r))]. (B.17) 
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In the following we show how to bring this equation into the form (1171) . We define 
a spectral function S{uj) fl20|) as the Laplace transform of the correlation function 
{5az{T)5az) ss of the qubit operators at the steady state. The steady state of the qubit 
should in principle be calculated with respect to the shifted Liouvillian £q, which is 
written using the value for a (1B.9I) . 

2ig\^jJ^)^ 

m ' 'n 

The second term is included into the detuning A = A + 45f^u;m/(ci;^ + 7^)(cr2)ss, but the 
correction is small, [gjuj^'^ <^ 1. Therefore it is justified to approximate the detuning 
by A ?^ A and to calculate the steady state with respect to C^. The master equation 
for the motion of the resonator finally gets the form, 

p^ = (£- + £-)p^ 

- 5f^ [S{uj^) [a + a\ [a, pj] + S{-uj^) [a + a'^ , [a^ , pj]] 

- g'^{S{uJrn) - S*{-uj^))[a + a'^,p^a\ 

- g\S{-uJ - S*{uj^))[a + a\ p^a^]. (B.19) 

After applying the rotating wave approximation, where we neglect the fast oscillating 
terms which are an order g'^/{uj^J?) smaller than the slow terms, we obtain, 

p^ = -iUraWa,prn] - ig^ [im {^(c^m)} [o^a, Pm] + Im {S{-uj^)} [aa\ p^a]] 

+ {g^Re{SicOm)} + 7m(iV^ + 1)) I?H(Pm) 

+ (/Re{5(-^^)} + 7miVm) V[a^]iPm), (B.20) 

which can also be written on the form flT7|) by introducing the effective frequency uj^ 
( TT8l) and cooling and heating rates A± (iT9l) . Note that we discuss the results in section 
m in the frame where the resonator is displaced by the coherent shift a (IB.9P . The 



master equation in the unshifted frame contains additional terms ~ aa^,a*a which do 
not contribute to the cooling equation, and a term lap which shifts the final occupation 
number. 

Appendix C. Spectrum 

In this section we first derive a general expression for the force spectrum 5'(a;) given 
in equation (!20!) . in terms of the parameters F, Fd, A and fi of the qubit dynamics. 
Further we evaluate the real part of the force spectrum in the limit of well-resolved 
peaks, corresponding to the resolved sideband limit for ion cooling [33]. To begin, let 
us write the qubit dynamics flT^ in terms of the corresponding Bloch equation 

dt{a) = A{a)-f (C.l) 

for the Pauli operators {ax,o'y,az), with 

-T-Td A 0\ [ ^ \ 

A= \ -A -T -Yd -n , f = 

n -2F J V 2r y 



(C.2) 
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and steady state solution (ct)ss = A~^T . We further note that the spectrum in equation 
f l20|) is given by the Laplace transform of the correlation function (5cr^(t)(5cr^(0))ss for 
the fluctuation operators 6(t = a — (fT)ss, 

S{u) = {6a,{s = -tu)6aM)ss, (C.3) 

which can be derived from the equation of motion for {Scrz{t)6crz{0))ss- In a first step 
we derive the equation of motion for {6a) from the Bloch equations ( JC.ll) . 

dt{Sa) = A{6a). (C.4) 

Then using the quantum regression theorem [46] , we find the solution for the spectrum, 

S{Lo) = -{0,0,l)[iujl + A]-^B, (C.5) 

with 

/ -i{ay)ss - {a^)ss{o-z)ss 
B = {6a6a^)ss= i{<yx)ss- {o-y)ss{<^z)ss ■ (C.6) 

In principle, we have now arrived at an expression for the force spectrum in terms 
of the qubit dynamics through the matrix A and the vector B. Let us proceed one 
step further and write the inverse matrix in the spectrum (JC.51) as [iul + A\~^ = 
Adj[ico'l + y4]/Det[ia'l + A], where Adj(a;) denotes the adjugate matrix of x and Det(x) 
denotes its determinant. The spectrum then gets the form, 

S{uj) = /"[""^ h{uj) = -(0,0,l)Adj[z^l + A]5, (C.7) 

where e^ are the eigenvalues of the matrix A. From this last step it is clear that the 
spectrum generally exhibits three peaks, with position and width given respectively by 
the imaginary and real parts of the eigenvalues of A. 

For our purposes it is sufficient to look at the real part of the spectrum ^(ci;), as 
the imaginary part only contributes a small shift to cjm, but does not contribute to the 
cooling. In the following we focus on the case of well-resolved peaks, which intuitively 
assumes weak qubit decay F, T^. With a spectrum containing only a few narrow 
peaks, cooling is significant only when the peaks are situated at the trap frequency 
cUm, and we therefore evaluate the spectrum only in their vicinity. To lowest order in 
F/A, Fd/A <^ 1, with A = v^M-A^, the e^ are approximated by, 

eo = -Fo, e± = ±tA - F+, (C.8) 

with effective decay rates Fq, F+ given by. 



and 



K^^;j^)(-H^:^"- 
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Thus in this hmit there are three well- resolved peaks at cu = 0, ±A. Cooling is optimized 
for A = tUm and close to the two peaks at ±A the real part of the spectrum can be 
written in the form ( 12^ with a± generally given by, 

"--^{ |±.A-r!|gk-r,| }- <«•"' 



which to lowest order in F, Fd <^ A are described as in flMj) . Note that the amplitudes 
of the peaks are given by a±/Tj^. The analytical expression for the real part of the 
spectrum will be evaluated in section H] with focus on the limit Fj/F -^ where we can 
neglect the contribution of the pure dephasing to the qubit dynamics. For finite pure 
qubit dephasing the cooling rate W (1261) generalizes to 



~ ~ ^siii (/JA/ I — Sir 

ly = 27, 



4 sin^ V^ V 1 — sin^ Lp 



J (4 - sm^ Lp) + (Fd/F) (2 - sm^ ^p) 

with siny? = Vt/K. The corresponding finite occupation number {n)^ (l30l) reads 
,, A^„|«»+(.p) + (rj/r)tt>_(.p)|+,3a_ 2±smV 



with a± given in ( 1241) . 
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